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Abstract 

We develop a method to analyze the strong coupling limit of the Bethe ansatz equations 
supposed to give the spectrum of anomalous dimensions of the planar Af = 4 gauge 
theory. This method is particularly adapted for the three rank-one sectors, su(2), su(l\l) 
and si (2). We use the elliptic parametrization of the Bethe ansatz variables, which 
degenerates to a hyperbolic one in the strong coupling limit. We analyze the equations for 
the highest excited states in the su{2) and sm(1|1) sectors and for the state corresponding 
to the twist-two operator in the sl(2) sector, both without and with the dressing kernel. 
In some cases we were able to give analytic expressions for the leading order magnon 
densities. Our method reproduces all existing analytical and numerical results for these 
states at the leading order. 
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1 Introduction 



There is mounting evidence that the planar M = 4 SYM gauge theory is integrable 
[1, 2, 3] to higher orders in perturbation theory, nourishing the hope to prove by this 
mean its equivalence with the string theory on AdS^ x S 5 [4, 5, 6]. Under the assumption 
of integrability, the spectrum of the anomalous dimensions of the planar Af = 4 SYM 
gauge theory is encoded, at least for states with large number of operators, in the Bethe 
ansatz equations. The S'-matrix behind these equations [7] is fixed by symmetry [8, 9], 
up to a scalar factor which is supposed to obey a crossing symmetry relation [10]. At 
strong coupling, the scalar factor should allow to reproduce the string theory results. Its 
first two orders in the strong coupling expansion were fixed in this way [11, 12]. More 
recently, based on previous work of Beisert, Hernandez and Lopez [13], Beisert, Eden and 
Staudacher [14] proposed an expression for the scalar factor which is crossing symmetric 
and reproduces both the strong coupling, string results and the perturbative results for 
the anomalous dimensions in the gauge theory. An efficient way of testing the predictions 
of the Bethe ansatz equations proved to be via the anomalous dimension of the twist-two 
operator. In the regime of large Lorentz spin S, this quantity scales logarithmically 

A-S = f(g)lnS + ... . (1.1) 

The universal scaling function for A/" = 4 SYM up to third loop order was extracted [15] 
from a perturbative three-loop computation in QCD [16]. The same quantity appears in 
the iterative structure [17] of the multi-gluon amplitudes of the supersymmetric gauge 
theory. It was computed up to three loop order in [18], and, by an impressing effort, up 
to four loop order in [19, 20]. The four loop result agrees with the prediction based on 
the proposal [14]. On the other hand, at strong coupling the universal scaling function 
f(g) is predicted to behave as [21, 22] 

f(g) = 4g- — + ... ■ (1.2) 

7T 

Here we adopt the recent convention for the coupling constant g 

? = ^ ■ (1.3) 

After the proposal [14] was made, several groups [23, 24, 25] attempted to derive from 
it the strong coupling expression (1.2). While the numerical work [24] easily reproduced 
the behavior (1.2) and even predicted the next term in the expansion, the analytical 
treatment proved to be much more difficult. Up to now, only the first term was obtained 
analytically [23, 25]. 

One of the aspects which deserve further attention is to understand the origin of 
the dressing phase [14] and to test its validity. It is believed that the dressing phase 
comes from a non-trivial structure of the vacuum. Very recently, a structure similar to 
that of the dressing phase was obtained [26], via the nested Bethe ansatz, for one of the 
non-trivial "vacuum" states in one of the sectors which are not of rank one. 

The aim of the present paper is to develop a systematic approach for solving the 
Bethe ansatz equations at strong coupling. Our analysis is based on the formulation of 
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the Bethe equations as integral equations, and it works most simply on, although it is 
not restricted to, states as the antiferromagnetic state in the su(2) sector, on the most 
excited state in the su(l|l) sector, or for the finite-twist operators, in the sl(2) sector. 
Some of this cases were already studied numerically or analytically; here we aim to a 
unitary treatment and hope that this method will be useful to obtain systematical 1/g 
expansion. 

Our starting point is the elliptic parametrization of the variables u and x ± that 
appear in the Bethe equations. The elliptic parametrization we use here is related to 
the one proposed by Janik [10] by a Gauss-Landen transformation, and it appeared 
independently in [9] and [27]. The elliptic modulus is defined by 



In the hyperbolic limit g — > oo, this parametrization supplies natural variables in the 
different regimes, recently discussed in [13, 28]. 

As it became clear from the numerical solutions of the integral equations [29], the 
points u = ±2g become singular for large g and they split the real axis u into two regions, 
with \u\/2g less or larger than one, where two different hyperbolic parametrizations 
apply. 1 The region \u\ < 2g is that of the "giant magnons", in the language of [28], 
and is characterized by finite values of the periodic momenta p. The region \u\ > 2g 
corresponds to the "plane- wave" limit [32], where the momenta are of the order p ~ 1/g. 

It is more difficult to see what happens exactly at the points u = ±2g and in their 
vicinity. These points separate the two regimes discussed above, and they are not prop- 
erly described by either of the two parametrizations. In fact, they correspond to momenta 
of the order p ~ g^ 1 ^ 2 , or to the "near-flat space" region, again in the terminology of [28]. 
In some cases almost all the roots of the Bethe equations are concentrated in this region. 
It is possible to use the elliptic parametrization to obtain the strong coupling expansion 
in this regime, 2 and it is clear that the expansion in this region involves powers of g^ 1 ^ 2 , 
or equivalently 1/A 1 ^ 4 . 

In order to obtain the strong coupling limit of the dressing kernel [14] we have several 
options. One is to take the large coupling limit term by term in the series defining the 
dressing phase around g = oo [13, 14]. For the giant magnon and plane- wave limits this 
works, and this is probably the most straightforward way to obtained the leading orders. 
In the near-flat space limit, however, as it was noticed in [28], all the terms of the series 
contribute to the leading order in 1/g. We would like to emphasize that, although the 
results for the leading order of the anomalous dimensions can be obtained without being 
particularly careful about the near-plane wave regime, the corrections will be crucially 
determined by it. 

A second possibility is to use the representation of the dressing kernel given by the 
"magic formula" of [14], and this representation was already exploited in [25]. In this 

1 The two hyperbolic parametrizations are related by the transformation s — > K — s, where K is the 
real period of the elliptic function. 

2 In this case, the parametrization is obtained by shifting s — > ±K/2 + s. 
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paper we give another related integral representation 

/oo 
dv Kb n (u,v) {Kl>-\v,u') - K^\v,u')) , (1.5) 



where the kernels K±(u,u') are the odd/even parts of the basic sw(l|l) kernel K(u,u f ), 
the inverse Fourier transform counterparts of the kernels Ko^(t,t') in the BES decom- 
position [14]. The upper indices K±' n (u,u') mean essentially that the variables u and u' 
are defined by elliptic functions of modulus 

k'n _ n Kn _ Q q\ 

K 40 ' k m Ag' 

Similar variables were used in [13] and they are associated to bound states of magnons 
[35, 36]. The sum in (1.5) resembles to a sum over intermediate states; in particular, the 
energy associated to these "intermediate states" is of the type [35] 

E n = sjn 2 + 16# 2 sin 2 p/2 . (1.7) 

In the large g limit, the sum over n can be taken with the help of the Euler-Maclaurin 
formula and the integrals over the variables v and z = n/4g are simple enough to be done. 
After taking the integrals, we reproduce the first term in the strong coupling expansion 
of the dressing kernel [11]. Although technically this result is not very useful, it can be 
considered as another evidence that the weak coupling and strong coupling definitions 
of the dressing kernel [14] agree. 

Turning to the explicit results, we obtain here the leading order result for anomalous 
dimensions of the highest excited state in the su(2) and sw(l|l) sectors, both for the 
Bethe ansatz with the Beisert, Hernandez, Lopez/Beisert, Eden, Staudacher (BHL/BES) 
phase, and without it. Even if the kernel without the dressing phase is by now only of 
limited interest, we can still ask the question whether the su(l|l) sector arises from an 
effective model, in the same way the su(2) sector appears from the reduction of the 
Hubbard model at half-filling. The large g limit of the su(l|l) sector should give a hint 
for the underlying degrees of freedom. In the su(2) sector, which was solved in [37, 38], 
the degrees of freedom are magnons, which become free as g — > oo. Our result in the 
sw(l|l) sector shows that these degrees of freedom do not correspond to a free system, 
and they are rather difficult to characterize. Two thirds of the excitations are of the 
type "giant magnon" , with finite momenta, and one third are of the type "plane wave" , 
with momenta of order 1/g, see also the numerical results of [39]. We find for the energy 
without the dressing factor 

E su (i\i) = —^—gL , E su( 2) = —gL . (1.8) 

The su(2) energy is the g — > oo limit of the exact result [37, 38]. When taking into 
account the BHL/BES phase, almost all roots of the Bethe equations, both for the su{2) 
and sw(l|l) sector, fall into the region "near-flat space" regime, with momenta of order 
l/y/g. This phenomenon was first discussed in [11], where the A 1 / 4 behavior was obtained 
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from the BA equations. The integral equations are able to reproduce the leading result 
for the anomalous dimensions of the corresponding states, which were also obtained 
numerically and analytically in [39, 40] 

EL m = V^9L, E d su{2) = ^L. (1.9) 

We also analyzed the strong coupling limit of the equation derived by Eden and 
Staudacher (ES) [29] for the anomalous dimension of the twist-two operator without and 
with the dressing phase (BES equation). Similarly to the approaches [24, 23], we obtain 
that the strong coupling limit of the ES equation is pathological. The BES equation has 
a better behavior, and we are able to reproduce the leading term of the density obtained 
recently by Alday et al [25] using the Fourier space representation, 

a< ( u ) = ^ for \u\<2g, (1.10) 

cf > {u) = - - - (1 — cosh s/2) for w = 2g coth s , \u\ > 2g . 

This solution reproduces the leading behavior of the universal scaling function (1.2). 

For the next order in the density, an important fraction of the roots lie near the points 
u ~ ±2g. We show that without properly considering the scattering of the magnons in 
this region, the density is non-normalizable. We leave the fine analysis of the near-flat 
space region for future work. 

The paper is organized as follows: in section 2 we write down the integral equa- 
tions corresponding to the three sectors, in section 3 we present in detail the elliptic 
parametrization and we take the strong coupling limit of one of the the building blocks, 
the sw(l|l) kernel. In section 4 we compute the energy of the highest states for the su(2) 
and sw(l|l) without the dressing kernel and discuss the leading order solution of the 
Eden-Staudacher equation. In section 5 we give an integral representation of the dress- 
ing kernel and we check that it reproduces the already known results at strong coupling. 
We also give the strong coupling limit of the AFS phase [11]. The integral equations 
with the dressing phase equations are considered in section 6. 



2 Bethe ansatz equations in integral form 

If we consider highly excited states, with a large number of magnons, the Bethe ansatz 
equations can be solved by transforming them into integral equations. These integral 
equations are particularly simple when the number of magnons is maximal, such that 
there are no holes in the magnon distribution 3 . Further simplification arise if we restrict 
ourselves to one of the three rank-one sectors, sw(l|l), su(2) or sl(2). The su(l|l) 
comprises one complex boson Z and one fermion U, and the state with the maximum 
number of fermions 

TrU L 

3 The integral equation for the highest state in the su(l|l) sector was written down by one of us and 
M. Staudacher, [41] and in [44], the equation for the antiferromagnetic state in the su(2) sector was 
solved in [37], and the equation in sl(2) sector, or the Eden-Staudacher equation, was derived in [29]. 
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is, at least at weak coupling, the one with highest energy. 

In the su(2) sector, corresponding to two complex bosons Z and X, the state with 
highest energy is, again at small coupling g, the anti-ferromagnetic state. This is, in the 
case of even length, the singlet state 

Tr(ZX) L/2 + ... . 

The sl(2) sector corresponds to combinations of one boson Z and the covariant derivative 
D. The interesting operators are the so-called twist-L operators with spin S, 

Tr D S Z L . 

The equations in the three rank-one sectors can be written in the compact form [8, 30] 

4-\ L if , 2 / + — / + —\V 

rr* 1 \ I {~t I If* If* I If* If* \ 



where r\ = 1, 0, —1 for su(2), sw(l|l) and sZ(2) respectively. The variables x ± are defined 
by ^ 

i« ^1 + v 7 ! - V/m 2 ) , x ± = x(u ± i/2) . 



X[U) 



It is often convenient to work with the momentum p, the physical variable of the magnon, 
defined by 



u{p) = 2 C0t |y 1 + 16^ 2 sin 2 | , or x + /x~ = e ip . (2.2) 

We are going to consider both the case where the phase 8(u k , u{) is zero and the case 
where it equals the BHL/BES dressing phase [13, 14]. 

Generically, the Bethe equations for the rank-one sectors can be written as 



M 



■Q gitfufc.u,) ? (2.3) 



e 

i=i 

where <p(itk, u{) is the scattering phase for two magnons. Let us exemplify the derivation 
of the integral equation on the su(l|l) case, which is the simplest one. We consider for 
convenience L odd. The total number of magnons is M = L. When g — 0, the equations 
(2.1) correspond to a free fermion system, with the occupation numbers 

P <«»=2rf/L, k = -±Z±,...,±Z±. (2.4) 

When g increases, the momenta start to evolve, according to the Bethe ansatz equations 
(2-1) ' 

|4 + i&^^)' ( 2 - 5 ) 



=i 
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When the length of the chain is very large, we can take the continuum limit. We introduce 
the variable t — k/L and the density of rapidities p{u) = —dt/du. The derivative with 
respect to u of the equation (2.5) gives 

i dv r°° 

p{u) = - — -£ + J dv! K(u, u') p(u') , (2.6) 
where the kernel is the derivative of the scattering phase 

The derivation of the su(2) equation is similar, with the difference that the maximum 
number of magnons is L/2 instead of L. This derivation follows closely the analysis 
of the antiferromagnetic state in the XXX model [42]. A more sophisticated treatment 
which is able to take into account the finite size corrections, was used in [43] to study 
the antiferromagnetic state in the BDS model [30] and in the Hubbard model. 

The Bethe integral equations (2.1) can be formulated in terms of three basic kernels. 
The first one is the kernel K(u, u'), corresponding to the su(l|l) case without the dressing 
phase. The second one is the su(2) kernel K su (2)(u,u'), and the third one is the dressing 
kernel, corresponding to the dressing phase 0(u,u') 

K(u,u') = --fflnfl- ^ f , 7 t-r, X \ (2-8) 

V ' ; 2mdu\ V x+(u)x-(u')J V x-(u)x+(u')J J V ; 

. 1 d , / u — u' + i\ 1 1 , , 

= ^-l-ln - - t =-- T- ., A2 , = (2.9) 



2tx i du \u — u' — i J n (u — u') 2 + 1 
K d (u,u') = L±e{u,u'). (2.10) 

When the dressing phase is taken into account, the complete kernel is 

/Ctot («,«') = (1 -v) K (u,u') +V K su(2)(u,u') + K d (u,u') 

= K,(u,u') +r}K su ( 2 ) («,«') ■ (2.11) 

It will be sometimes convenient to work with the phase associated with the kernel /C 

K(u, u') = ^du^u, u'), = (1 - v)<P + 0. (2.12) 

The sl(2) sector is special in the sense that the number of magnons, also called spin, 
M — S, is not bounded, even if the length is finite (L = 2 for the twist-two operator). 
The lowest state with S magnons can be characterized by the numbers 

n k = k + sign(/c) , k = -(S - l)/2, . . . , (S - l)/2 . (2.13) 

The corresponding integral equation is given, again in the absence of the dressing phase, 
by 

p(u) = -5(u) + - / du' P{ + 2 / du' K(u, u') p{u') , (2.14) 

S 7T (U - U'Y + 1 
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where the term proportional to 5(u) comes from the distribution of the mode numbers 
nfc. The limit g = of this equation is singular, and it was solved in [31]. The solution 
for large S is 

1 , 1 + a/1 - Au 2 /S 2 2\nS 

p " {u) --s^ i-^-J /S ^ — + ° {llS) - (2 - 15) 

Eden and Staudacher [29] chose to separate the density p(u) into the g = part, po(u), 
and a fluctuation a(u) 

p{u) = Po ( u )-2g 2 ^- a(u)+0(l/S) = -2g 2 ^- ( a(u) - +0(1/S) . (2.16) 



Aug 

From (2.14), the equation satisfied by the density fluctuation a(u) is 

«*" = i jy j^vh + 2 jy k{u ' u,) h - *y • <2 - i7) 

This is essentially the equation derived by Eden and Staudacher [29], with the inhomo- 
geneous term written differently. This way of writing has the advantage to show that the 
separation of p(u) into po(u) and a(u) makes sense, at large g, if K{u,u') is sufficiently 
well behaved at large u' . As we shall see, if K(u,u') = 0(1) for large g and \u'\ > 2g, 
the equation (2.17) will have only the trivial solution 4 , since at leading order in 1/g one 
can write 

J du' K(u,u') (a(u') - = . (2.18) 

In this case, the separation of the density into a free part and a perturbation is inconsis- 
tent. The Eden-Staudacher kernel is exactly of this type, and this might be the reason 
the strong coupling limit of the Eden-Staudacher equation is so pathological [23, 24]. 
The BES kernel, on the other hand, vanishes at the order 0(1) for \u'\ > 2g, and the 
strong coupling limit of the solution of the BES equation is well behaved. 



3 Elliptic parametrization of the kernels 
3.1 The elliptic map 

In the strong coupling regime, it is convenient to rescale the variables u and x ± in order 
to eliminate the overall factors of g. We will use the rescaled variables throughout the 
rest of the paper 

U = ^, X ± = ^, €=—. (3.1) 

The function 

x{u) = u + Vu 2 - 1 (3.2) 

has a branch cut along the interval [—1,1]. Since the Bethe equations involve the shifted 
variables u ± ie, they contain two symmetric cuts, [— 1 — ie, 1 — ie] and [—1 +ie,l + ie\. 



'At least if the kernel is non-degenerate. 
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The cuts can be removed by introducing a global parametrization. We eliminate the 
parameter u from 

u ± = u±ie = \(x ± + l/x ± ) (3.3) 
to obtain the following relation between x + and x~ 




If we require that the momentum p is real, then x + and x~ are complex conjugate and u 
real. Then the condition (3.4) for x = x + and x = x~ define a contour C in the complex 
x-plane. The contour consists, for g real, of two connected components, and Co, 
which are exchanged by the particle-antiparticle transformation x — > 1/x. We denote 
by Coo the component that contains the point x = oo; the other component contains the 
point x — 0. In the limit e — > the contours and Cq develop cusps where they touch 
(Fig.l). At this point the contour C is the union of the real line and the unit circle. A 
second singular value is e = ±i. Here the two connected component join into one (Fig. 2). 




e > e = e < 



Fig.l : The contour C in the x plane for three real values of e close to 0. 




ie < 1 ie = 1 ie > 1 



Fig. 2 : The contour C for three imaginary values of e close to —i. 

The equation for the contour C can be written conveniently in terms of the momentum 
p and a complementary parameter (5 

x ± = e (P±i P )/2 (35) 
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In these variables the relation (3.4) reads 



sin - sinh — = e (3.6) 

and the parameter u is expressed as 

u = cos ^ cosh ^- . (3.7) 

In the following, we consider only positive values of e, and p G [0,27r]. This insures 
that u in (3.7) takes all the real values once and only once. 

The relation (3.6) is resolved by elliptic map 5 with modulus k defined as 

k! 1 
e = — or k — . . (3.8) 

Then f3 and p are parametrized by the Jacobi elliptic amplitude function 

p(s) = 7T - 2am(K - s, k), @(s) — —in — 2i am(iK' — s,k). (3.9) 

When the elliptic parameter s sweeps the interval [0, 2K], the momentum p increases 
from to 27T. The symmetries of p and f3 are 

p(-s) = -p{s), p{s + IK) = p(s) + 2tt (3.10) 

P(-s)=P(s) + 2m, (3(2K-s) = (3(s). (3.11) 

The functions we will work with are actually expressed in terms of Jacobi elliptic func- 
tions 



e 



±/3/2 



1 ±dns 



, e ±ip/2 = cds± ik' sd s 



ksns 

d s p = 2k'nds, d s (3 = -2css. (3.12) 
The contours and C are parametrized respectively by s G 1^ and s G 1$, where 

Joo = [0, 2K], J = [2iK', 2K + 2iK'\. (3.13) 

The original rapidity variable u is parametrized as 

,, 1 ens 
k sn s dn s 

The map (3.14) obeys the symmetries u(s) = —u(—s) = u(s + 2K) = u(s + 2iK'), as well 
as 

u(s)u(K-s) = 1 + e 2 . (3.15) 



'Our map is related to that of [10] by a Gauss-Landen transformation fcj an ik = 2v / fc' /(l + fc')- 



The real axis in the w-plane is the image of the interval 1^. The functions x (s) are 
periodic in s with periods 2K and 4iK' and have the symmetries 

x +( s ) = x+(K-zK'-s) = -x-(-s), x ± (s + 2tK') = (3.16) 

x ± (s) 

The first symmetry preserves the contours and Co, while the second symmetry ex- 
changes them. Note that the factors e 13 ^ 2 and e ±ip ^ 2 are themselves antiperiodic in 
s — > s + 2K. This is a little complication, due to the fact that we are working with 
trigonometric functions of p/2, which are periodic with period Air. However, in the phys- 
ical variables, they always appear in combinations which are periodic with period 2tt. 
As we already mentioned, we prefer to work with p in the interval (0, 2tt). 




In the limit e — > 0, where K — ► oo and K 7 — ► n/2, the elliptic parametrization de- 
generates 6 and the variables x and u become hyperbolic functions of s. As the plot of 
u(s) in Fig. 3 suggests, the parametrization interval splits into two different pieces, cor- 
responding to two different regimes. By periodicity 7 , we can choose the parametrization 
interval as [— ^K, |K]; then the first regime corresponds to the interval [— |K, |K], which 
becomes the whole real axis when K — > oo. The second regime, corresponding to the 
interval [|K, |K] can be mapped also to the interval [— |K, |K] by the transformation 

s^K-s. (3.17) 

The first interval corresponds to values \u\ > 1, while the second corresponds to \u\ < 1. 
We therefore introduce the suggestive notation 

!> = [-§K,§K], I< = [|K,§K] (3.18) 

6 Another point where the torus degenerates into a cylinder is e — ±i, or g 2 — —1/16, cf. Fig. 2. 

7 As already explained, the variables u and x are periodic with period 2K, and we can safely do the 
shift s — > s — 2K for them, but for the variables p, /3 we have to remember that the point — iK was 
originally |K. 
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and use the asymptotic formulas valid in the limit e — > and for s G I> 

sns = tanhs, sn (K — s) =1 

ens = l/coshs, en (K — s) — k'sinhs, 

dns = l/coshs, dn (K — s) =k'coshs. (3.19) 

The next terms in the expansion in e are given in Appendix A. 

In the hyperbolic limit, the expression (3.14) for u becomes, in the two different 
regimes, 

■u > (s) = u(s) ~ coths , «<(«) = u(K — s) ~ tanhs . (3.20) 

The points s = ±K/2, which are the fixed points of the transformation s — > K — s, 
correspond to the points u = ±Vl + e 2 - In the strong coupling limit, e — > 0, these 
are the points -u = ±1, and they correspond to meeting points s — > ±oo of the two 
parametrizations (3.20). These points will be very important, and we need a better 
parametrization for them. In the strong coupling limit the vicinity of these points is 
parametrized by (Appendix A) 

u(±\K + s) = ±1 - e sinh 2s + C(e 2 ). (3.21) 

The expressions for the momentum p and its counterpart f3 in these regimes are 

sinh/3/2 = v^e Ts sinp/2 = ^~ee ±s , (3.22) 

where the upper sign correspond to momenta close to and the lower sign to momenta 
close to 2n. A similar parametrization was used in [28]. 



3.2 Strong coupling limit of the sw(l|l) kernel 

We saw that in the strong coupling limit the integration interval splits naturally into 
two domains, \u\ < 1 and \u\ > 1, which are the images of the intervals I< and J> in the 
s parametrization. As it is clear from Fig. 1, in the first domain the complex variable 
x = x + becomes asymptotically unimodular, while in the second domain it becomes 
asymptotically real. For e — > the relations (3.6) and (3.7) give 

-> , u = cosp/2 if \u\ < 1 , (3.23) 
p^O, u = cosh/5/2 if \u\ > 1 . (3.24) 

Passing to the elliptic parameter s we write the integral equation (2.6) as 

P( s ) = 7^P(s) + I d Sl K(s, Sl ) , (3.25) 

where 

p(s) = \u'(s)\p(u), K(s, Sl ) = \u'{s)\K(u, Ul ). (3.26) 

The first term in the r.h.s. of (3.25) changes sign because the derivative u'(s) is negative. 
The kernel is equal to the derivative 



K(s, Sl ) = - — 4- flnfl-e-^+^) e -^(p-PiA -lnfl- 
v ' ; 2irids I \ ) V 
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J_ l( ) _ 1 p'fc) smh |(g + gi) - (3'{s) sin |(p - Pl ) 

Let us mention here that the piece p'(s)/4ir has a physical meaning: it corresponds to a 
change of the effective length of the chain. Every magnon increases the effective length 
of the chain by 1/2. This is particularly clear on the discrete equations, where, at strong 
coupling, the Bethe equations become, as it was noticed in [44] 

e i Pk (L+M/2) = 1 ^ (3>28) 

under the condition that all p^ are finite. 

It is now easy to take the strong coupling limit using the asymptotic expressions for 
p(s) and j3(s), which can be found in Appendix A. The asymptotic form of the kernel 
is given by four different analytic expressions depending on whether its arguments are 
in the interval I> or I<. After applying the redefinition (3.17) to the arguments that 
belong to the second interval, we can write the result as a 2 x 2 matrix kernel 

K»(s, Sl ) K><(s, Sl )\ (K(s, Sl ) K(s,K- Sl ) \ . , 

K<>(s, Sl ) K«(s,s 1 )J • \K(K-s, Sl ) K(K- s,K- Sl )J ' 1 ' 

The arguments s,Si are defined in the interval I> , which extends to [—00,00] when 
e — > 0. Similarly the density splits into two components, 



pHs)\ = ( p{s) 
p<(s) U(K-a) 



(3.30) 



satisfying, in the limit e — >• 0, the normalization condition 8 

/oo roc 
dsp < + / dsp > = 1. (3.31) 
-OO J —CO 

Eq. (3.25) now takes the form 

/oo noo 
ds 1 K»(s,s 1 )p>(s 1 )+ / ds 1 K><(s,s 1 )p<{s 1 ), (3.32) 
-OO J —OO 

P < (s) = ^+ [°° ds 1 K<>{s,s 1 )p>(s 1 )+ [°° ds 1 K«{s,s 1 )p<(s 1 ). (3.33) 
^ J -00 J -00 

The matrix elements of the kernel (3.29), evaluated in Appendix A, are 

K«(s,s>) = ±^—-5(s-s'), 
2tt coshs 



111 1 



K<>(s, s') 

K><(s,s') 

2n cosh(s — s') 

K»(s,s") = 0. (3.34) 



2rr coshs 2tc cosh(s — s') 
1 1 



3 This is true only for the kernel without the dressing phase, at the leading order. 
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Note that in the hyperbolic limit all four integration kernels depend on the difference 
of the arguments. This is a considerable simplification, since the equation can be now 
solved using Fourier transform. 

In the third regime, which corresponds to the points around u — ±1 the momenta 
are parametrized by (3.22). This regime will be very important for evaluation of the 
density with the BHL/BES kernel. The appropriate parametrization is given by (3.21) 
and (3.22). The only component of the undressed kernel that survives at the leading 
order is when either u,u' ~ 1 or —1 

K ++ (s, s') = K—(-s, -s') (3.35) 

= _ J_ sinh( S - S Q + l + ^ 

47r cosh (s — s') cosh (s + s') — sinh (s + s') 

If only one of the arguments is in the intermediate region, the kernel is of order \fe. 

One might argue that there are many intermediate regimes, where p scales like an 
arbitrary power of e, 

p ~ e 1 " 7 , ~ e 7 7 G (0, 1) . (3.36) 

However, if 7 ^ 1/2, one of the variables p and (3 dominates the other, and these regimes 
are properly taken into account in the regions <. Only for 7 = 1/2 the quantities p 
and (3 are of the same order of magnitude and this is why we have to consider this case 
separately. 



4 Solving the integral equations in the strong cou- 
pling limit without the dressing kernel 

4.1 The su{2) case 

The su(2) kernel is the only one which is of difference form, therefore it can be exactly 
solved for any value of g [37, 38]. The energy of the antiferromagnetic state is identical to 
the energy of the ground state for the Hubbard model at half filling [46]. Here, we would 
like to comment on the strong coupling limit of this solution. The integral equation is 
particularly simple, since the kernel becomes simply the delta function 

K sum (u,u>) = -K™l {u _ u ,y + Ae 2 = -Kn-W). (4.1) 
The integral equation reduces then to 

The solution corresponds to L/2 roots distributed uniformly between p = and 2tt. In 
this case, the density is normalized to 1/2, because the total number of magnons is L/2. 
All roots of the Bethe equations are inside the interval \u\ < 1 and the energy is 

E su{2) = AgL [ dp I sinp/2| p(p) = . (4.3) 
Jo 71 



13 



Physically, the excitations are the (usual) magnons. The solution obtained above 
shows that, in the strong coupling limit, the magnons are free, except for the statistical 
repulsion (4.1). The situation is very much similar to that of the Haldane-Shastry model 
[47, 48], where the magnons are also free up to the statistical repulsion, translated into 
the rule that two magnons cannot occupy successive momenta. Let us remind that the 
scattering phase for particles with purely statistical interaction is [49] 

(Pip, P') = (A - 1) 7T sign(p - p') , (4.4) 

which is consistent with (4.1), the statistical parameter of the magnons being A = 2. 
The only difference with the Haldane-Shastry spin chain is that here the magnons have 
the dispersion relation 

E(p) = Ag\ sinp/2| , (4.5) 

which is typical for a finite-difference hamiltonian, while the Haldane-Shastry magnons 
have the dispersion relation 

E HS {p)=p{2ix-p) . (4.6) 

A natural candidate for a model with purely statistical interaction and with trigonometric 
dispersion relation is the Ruijenaars-Schneider model [50]. This suggests the existence 
of a spin model which would describe the spin sector of the half-filled Hubbard model 
for any value of g and which might be a multi-spin deformation of the Inozemtsev model 
[51,52]. 



4.2 The su(l\l) case 

The solution of the su(l|l) sector is considerably more involved, but its leading order 
still can be obtained in closed form. Since now the roots of the Bethe equations occur 
both in the regions \u\ < 1 and |u| > 1, the integral equation splits up into two coupled 
equations for the densities in the two regions 



2p<{s) 



2n cosh 



1 -'f^ (4 . 7) 
shs 2tc J^oc cosh(s — s') 



p>(s) = -L ds' f/^ r . (4.8) 
2tt 7_ 00 cosh(s - s') y ' 

Substituting the second equation into the first we obtain a new equation where the kernel 
is of difference form 

3 1 1 f°° s — s' 

P < ( s ) = ~A T~ - T"2 / ds ' -T7 7\ P K ( s ') > ( 4 - 9 ) 

47rcoshs 47r 2 7_ 00 smh(s - s') 

and therefore is solved by Fourier transformation. In Fourier space the equation becomes 
algebraic, 9 

3 7T 1 7T 2 

P<it) = 47cosh(vrt/2) " 4^2coshV/2) P<(t) ' ( ] 



9 In order to keep the notation simple, we are going to use the same symbol for the Fourier transforms 
p< (t), p> (t) of the densities. 
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with the solution 



Transforming back to the s-variable, we get 

/ x 1 cos2as/n . t . 

p< (s) = -= -1— , 4.12 

v/2tt coshs 

a = lnv^ ~ .34657 . (4.13) 

For the density of roots p{u) we obtain in the interval \u\ < 1 

/ x P < ( s ) 1 1 /'a, l + ,,,,, 
p< (u) = f-^L = - cos - In . 4.14 

' ^ \U'{S)\ ^ y/T^? \n 1-uJ V ^ 

To evaluate the function p > (u), we make the change of variable u = coths, which is 
appropriate for \u\ > 1, so that we get 

> 1 f°° ds' cos(2as'/7r) 1 sin(2as/7r) 

^ 2v / 2vr 2 7-oo cosn s> cosh(s — s') n sinh s 

or, in the initial variables, 

p>( u ) = I_^ =sin (^i I1 ?i±i > \ . (4.16) 

The result shows a singularity at the points u — ±1; moreover, it can take negative 
values. This happens at values of 1 — u ~ 10~ 5 and is presumably due to the roots around 
u — ±1 which were not taken properly into account. Except for these oscillations, the 
result reproduces nicely the solution obtained by numerical integration of the Fredholm 
equation 10 . 

The total normalization of the density is insured, since J du p<{u) = 2/3 and 
f du p>(u) = 1/3. So two third of the roots of the Bethe equations fall in the re- 
gion | it | < 1 and one third in the intervals \u\ > 1. This phenomenon is also present in 
the numerical solution of [39]. Finally, the energy of the state is given by 



E. 



AgL j du P< {u) Vl-u 2 (4.17) 

2y/2gL f°° , cos2as/n 16a r . 

f ds J— = gL . (4.18) 



su(l|l) 



vr J_ 00 cosh s vr 

This result is consistent with that of Beccaria and Del Debbio [39]. We obtain for their 
constant cl the value cl = 21n(2)/7r 2 ~ .14046. 



10 



We thank M. Staudacher for supplying us the solution obtained by numerical integration. 
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Fig. 4: Density profile predicted by the analytic solution (4.14), (4.16). Only the u > part 
is shown. The oscillations do not show yet at the chosen scale. 

4.3 The s/(2) case 

The dimension of the twist two operator was originally supposed to be determined by 
the solution of Eden-Staudacher equation [29] 



a(u) 



7T 



du' 



2e 



u - u') 2 + 4e 2 



a{u') + 2 / dv! K{u,u') \a{u') 



2ng 



(4.19) 



The anomalous dimension is simply given by the normalization of the density fluctuation 
a{u) [53] 

/oo 
du a{u) , (4.20) 
-oo 

or alternatively by the formula [29] 

/(s) = tf (i . (4.2i) 

As explained in section 2, in the limit e = l/4g — > the equation (4.19) becomes 

1 



du' K(u, u') [ a(v!) 



2ng 



. 



Obviously, this equation has as solution 



ii 



<J(U) 



2ixg 



(4.22) 



(4.23) 



'Here we use the variable u rescaled by 2g, whence the extra factors of 2g in the density. 
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on the whole real axis. This solution, first obtained in [29], seems to be correct in the 
interval \u\ < 1, since it insures the vanishing of the 0(g 2 ) term in f(g) 

f(g) = Sg 2 (l - Ag du o< (u) + 0(l/<7)) = O(g) . (4.24) 

In the interval \u\ > 1, although it has the right scaling in g, the solution (4.23) is 
not integrable and it gives formally an infinite value for the anomalous dimension of the 
twist-two operator. As mentioned in the section 2, the density (4.23) exactly compensate 
the leading order in S in po(u) everywhere in u 

. , 4oln S , , ^ . 

A) « = — + - • ( 4 -25) 

7T O 

Unfortunately, as we show below, (4.23) is the unique solution to the ES equation at 
strong coupling. We can write the equations as 

J_ f ds > ^ = 1 I (4 26) 

2tt J cosh(s-s') 87r^cosh 2 (s/2) ' v ' ; 

1 f J , a > (s > ) <t . 11 1 A 

7T- Ids' / ^ - +<J < (s) = r^T + 



2n J cosh(s — s') Aug cosh 2 (s) 27rcosh(s) 

where A = f ds (o" < (s) + o _> (s)) is the normalization of the density. In Fourier trans- 
formed form, the equations become 

a<{t) t 



2 cosh nt/2 2g sinh nt ' 

ff> W +„<(«) * ■ A 



(4.27) 



2cosh7rt/2 45(sinh7rt/2 2cosh7rt/2' 
The first line is readily solved as 

1 11 1 

a<(t) = — — , or <t<(s) = 5- , or a < (u) = — , (4.28) 

w 2^sinhvrt/2' w 2ng cosh 2 s 2irg ' v ' 

while the solution for the second is 

a>(t) = A-- f - -, or a>(s)=AS(s)+ ' ! 



2gnt&nhnt/2 ' 27rg sinh 2 s ' 

or a > (u) = -A 5{u - oo) + — . (4.29) 

2ng 

The normalization A is not fixed by the equation. 

All other approaches tried to compute the strong coupling limit from the Eden- 
Staudacher ran into some pathology as well: the result obtained by Kotikov and Lipatov 
[23] oscillates strongly, the result obtained numerically in [24] may not converge to a 
straight line in g as expected 12 and the numerical solution for u{u) in [29] does not seem 
to converge to a value with definite normalization when g — > oo. We believe that all these 
pathologies are related to the separation of the density p{u) into po(u) and a(u). Rather 
surprisingly, this pathology disappear when the dressing kernel is taken into account, as 
we will show in the section 6. 



2 We thank the authors of [24] for checking the behavior of their ES solution. 
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5 The dressing kernel 



5.1 Integral representation 

Beisert, Eden and Staudacher proposed an expression for the dressing factor in the Bethe 
ansatz equations. This dressing factor translates into a correction to the kernel in the 
integral equation. Their formula is given in Fourier transformed form; for our purposes 
an expression in terms of variables u is more suited. Let us start with the Fourier 
transform of the su(l|l) kernel 13 

K(t,t) = -tt(1 - sign tt')\t\e-W + ^ J> . 

n>0 ' ' 

which is related to the kernel of [14] by 

K(t,t') = -|(1 - signtt')\t\e-^ + ^K m (\t\, \t'\) (5.1) 

and can be broken into a symmetric and antisymmetric part 

K(t,t') = K+(t,t') + K-(t,t') 

= signttOI^" (l ' l+l '' l)£ (^o(|tUt1) + ^i(|tUt'|)). 

In the ES and BES equations, the kernel is written in the Fourier variables defined 
only for positive values of t and t' . This is made possible since the densities p(u) and 
<j{u) are symmetric under u — > — u, and therefore under t — > —t, for the states under 
consideration. For generic states, this symmetry is not present, and this is why it is 
desirable to have a non symmetric kernel at hand. The extension of K d (t,t') to negative 
values of t and t' is ambiguous, and undoing the Fourier transform in the "magic" formula 
of [14] will be ambiguous as well. 

The authors of [25] gave a prescription to to define Kd(u,u f ) by symmetrizing all 
the kernels on the second variable K(u,u') — > ^(K(u,u f ) + K(u, —u')). We try to avoid 
symmetrizing the kernels from the beginning, and give here a different prescription, 
using instead that the full BES kernel is antisymmetric with respect to u± — > u T on 
each variable separately. We remind that the transformations u — > —u and u± —>■ u T 
translate in the following way on the variable p 

u — > -u ^ p — > 2ir -p (5.2) 

u± — > u T <^ p — > — p (5.3) 

while the variable (3 stays unchanged. Of course, these two transformations act differently 
on the Fourier transform. 

It will be simpler to work with the non-symmetrized kernels and impose the antisym- 
metry under (5.3) at the end of the computation. We define, for t,t' > 0, 

K d (t,-f) = -A7rte-^K c (t,t') 
13 Our definition of the kernel differs by a minus sign from the one of [29]. 
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poo ill 

= -^te'^j dt''k 1 ^t'')^ r — [ k {t\t') 

= J dt"K_(t, -n^p^K^", ~t') • 



For any values of t, t', we can define the dressing kernel as a convolution 

^(M') = — / dt"K_(t,t") 2 , n _ K + (-t",t') . (5.4) 

The dressing kernel K d {t,t'), defined in this way, is proportional to 1 — sign tt'. Note, 
however, that after anti-symmetrization under (5.3) this property will not hold anymore. 
Going back to the u variables we obtain 

/OO A pOO 

dv K-(u,v)K+(v,u') - - / dv dv'K_(u,v)h( y v -v')Q+(v',u') 
-co ^ J —oo 

(5.5) 

where 

, . . 2e f°° , \t\e itu 2e / 1 (vr/2e) 2 \ 

h(u) = — dt -fi = — I — V — " (5.6) 

KJ 27r7-oo e 2 ^l-l 2tt \m 2 sinh 2 (7TM/2e); 

and + (w,w') is the Fourier transform of K + (t,t')/\t\. The function h(u) becomes 5(u) 
as e approaches 0. Since the kernels involved in the convolution in (5.5) are of order e°, 
the dressing kernel will be of order e -1 . This formula is very close to the one in [25], 
except that we traded the difficulty in evaluating the function h(u) for the difficulty in 
evaluating the phase Q + (u,v). 

Let us now write down the expressions of K + (u, u') and K-(u, u'), which contain the 
odd, respectively even powers in the expansion of the logarithm 

K-(u,v) = - — :C> M ln 



with 



Am (1 -X_ + )(l + X_ + ) 

i i 

X+- = ~r\~rr > X-+ = . (5.9) 

x + (u)x (v) x (u)x + (v) 

A strategy to compute the expansion of K d in the strong coupling limit is to compute 

the limit of the auxiliary kernels (5.8), and then perform the integrals, as it was done in 

[25]. This procedure becomes rapidly quite involved, since different pieces of the kernel 

contribute for different orders in e. 

In the following, we give yet a different representation for the dressing kernel, which 

may be more useful for the strong coupling expansion. The equation (5.4) can be written 

as 

A . poo 

K d (t,t') = --Y dt"k_{t,t")e- 2n ^k + {-t",t'). (5.10) 
71 t?o J -°° 
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This formula can be easily Fourier transformed back such that we get 

/oo 
dvKb n (u,v) Kl'\v,u') (5.11) 



with 



K? n (u,v) = -±-d u ]((x [+m Hu)x(- n \v))- l -(x(- m \u)x( +n \v))- 1 ) , 

I~>0- evon 



X 



(±n), 



u) = u ± ine + y/ (u ± me) 2 — 1 . (5-12) 



The variables x( ±n ) appeared in [13], where they play an important role in defining the 
dressing phase. They live on a torus with modulus defined by 

K 

and they are associated with the bound states of magnons [35]. 

The expression (5.11) for the dressing kernel involves a sum of integrals along an 
infinite set of contours C™, n = 1, 2, 3, ... , depicted in Fig. 5, defined in the same way 
as the contour C°°, but with e replaced by ne 

= {x E C I (x - x) (1 - l/xx) = Aine, xx > 1} . (5.13) 




Fig. 5: The sequence of contours in the definition of the dressing kernel 
drown in the x-plane (left) and in the (p + i/3)-plane (right). 



In the strong coupling limit, the expression (5.11) becomes tractable, since we can 
evaluate the sum using the Euler-Maclaurin formula 

£/K) = i/ f^dz-lm-^-^f^Ho). (5.14) 

n=l €J ° 1 k=l [ZK) - 
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We thus extend ne to the continuous variable 

ne = z (5.15) 

and define the contour by 

Cf : sin | sinh | = z, (3 > 0. (5.16) 

The integrals over t> and the variable z become a double integral over the domain T> 
in the complex x-plane spanned by the contours C^°. The domain in question is the 
the upper half plane minus the unit semi-disk. In terms of the complexified momentum 
(3 + ip = 2 In x, the domain T> is the semi-infinite cylinder T> = {0 < p < 2ir, (3 > 0}. 

The integral can be done without specifying a particular parametrization of the con- 
tour Cf if we express the integrand as a wedge product of differential forms. We first 
write the basic kernel as a 1-form 

dp 1 sinh &&dp- sin ggjg 

A I.Ii =— h- ^-5— -5 . 5.17 

1 ' ; 4tt 4tt cosh cos » 1 ' 

Then the "even" kernel K + is a one-form in the basis of dp and 

1 sinh ^ cos ™ dp - sin ™ cosh d(3 

K + (x,x 2 ) — , r t r , (5.18) 

27r cosh(/3 + (3 2 ) - cos(p - p 2 ) 

while the "odd" kernel is a zero-form in this basis, 

_ 1 ^ dpxsinh^ + ^-^isin^-p) 

4vr^ 1 + 47r cosh^ + ^O-cos^-p) ' 15 9j 

Also, for the differential dz we have from (5.16) 

dz=\ sinh(/3/2) cos(p/2) dp + § sin(p/2) cosh(/3/2) dp . (5.20) 

Then the leading term in the strong coupling limit e — ■> is given by 

^(xi,x 2 ) = — dz K_K + = — / K„dzAK + . (5.21) 

The evaluation of this integral (the calculation is sketched in Appendix B) leads to the 
following expression for 

K° d (x 1 ,x 2 ) = —d u [-x(xi,xt)+x(xt,x];)+x(xT,-xt)-x(-x 2 ,x];)+c.c.], 

X (x,y) = -(^x + ^j\og(y-^j . (5.22) 

This expression must be antisymmetrized with respect to xf — > x^ (or, equivalently, 
p 2 — > — p 2 , (3 2 — > Z^)- After antisymmetrization the result coincides 14 with the AFS 
kernel K (u, v), symmetrized with respect to v — > —v. 



It is useful to notice that we can substitute xo with x to calculate AFS term in (5.23). 
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5.2 The strong coupling limit from the AFS kernel 

To compute the leading order(s) of the dressing kernel, we do not really need the integral 
representation derived in the previous subsection. It is more straightforward to use the 
strong coupling representation of the dressing phase 6\ 2 = 9(ui,U2) worked out in [13, 14]. 
The latter can be cast in the form 



1 



012 = ^ - Xt 2 ~ - Xu + + xt 2 + - (1 - 2)] , (5.23) 

Xu = X x s 2 ) , r,s = ±, 

where the function x can be expanded in powers of the inverse coupling constant 

X = Y J Xn{2e) n . (5.24) 

n>0 

To compute the leading term in the anomalous dimension, the AFS term [11] is sufficient, 

Xo(x,y) = ~ — In — , (5.25) 



V xy 

while for the next correction we will need the Hernandez-Lopez [12, 54] term 



Xi(x,y) 



7Y 



y+l x-y 



(5.26) 



+Li2 v^Vv^ _ u Vy + VVy + Li2 y/y-yyy _ Li2 Vv + 1 /Vv' 



^/y + Vx 



In the following, we give the first non-zero order of the dressing AFS phase. The kernel 
can be simply obtained by taking the derivative of the phase divided by 2n, with the 
right sign. We find, for the sectors " < " and " > " , 



0»(u 1 ,u 2 ) 

e><( Ul ,u 2 ) 
e<>( Ul ,u2) 
e«( Ul ,u 2 ) 



2e sinh (ft -p 2 )/4 

sinh ft/2 sinh fo/2 sinh(ft + p 2 ) /4 
2 sinp 2 /2 



-2^»( Ul ,u 2 ) + 0(e 2 ), 



sinh ft /2 
2 sinpx/2 
sinh ft /2 



-2<p><( Ul ,u 2 ) + 0(e), 
-2<p<>( Ul ,u 2 ) + 0(6), 



;cos Pl /2-cos P2 /2)ln S ^ (Pl - p2)/4 



0(1). 



(5.27) 



sm 2 ( Pl +p 2 )/4. 

We remind thet tp is the s«(l|l) phase. 

In the sectors "+" and "— ", where f3 and p (or 2n —p) are of order y/e, we introduce 
the parameter a, which is of order one, as follows 

Pi = 2\feai , ft = 2y/ea~ 1 with a« > in the region" + " , 
Pi = 2tt + 2\/eai, Pi = — 2y / ea i ~ 1 with ctj < in the region" — ". 
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Then we have 



6 ±T (ai, a 2 ) 
9 ±> {a 1 ,u 2 ) 

6 ±< {a 1 ,u 2 ) 
^(ui.aa) 



((Xi,Ct2) = ±77 {oil — «i — « 2 + Q 2 J In _ 1 



(a x 1 + a 2 x ) 2 + («i + a 2 ) 



-Ai In 



a l 1 + a 2 1 + ^(^1 ~~ ^2) 

±2 + C(v^) 



+ 0{V~e) , 



±2v^ 



sinh f3 2 /2 



+ 0{e) , 



■— ai sinp 2 /2 + C(l), 



T2v^ - 



a 2 



sinh ^i/2 



0(e), 



a 2 sinpi/2 + (9(l). 



(5.28) 



In (5.27) and (5.28) we omitted the signs coming from u being in u > 1 and it < — 1. 
A mnemonic rule to reproduce these signs is to take sign/3j = sign-Uj. In the hyperbolic 
parametrization with the variable s, which we give below, the signs are automatically 
taken into account 



s — s 



6»{si,s 2 ) = 2e sinh s sinh s' tanh ' L -^~ ~ 2ip»(s l ,s 2 ) + £>(e 2 ) 



0><{si,s 2 ) 
<> {s 1 ,s 2 ) 



2 sinh si 



cosh s 2 
2 sinh s 2 



-2(p><( Sl ,s 2 ) + 0(e) 



cosh si 



-2<f<>{si,s 2 ) + 0{e) 



a«< \ 2 / sinh(si - s 2 ) s x - s 2 , 

^ << (si,S2) = - — r 7 — In I tanh — - — 1+0(1) 

e V cosh s\ cosh s 2 2 



(5.29) 



In the regimes where one of the variables is in the region " + " or " — " , we put 

ai = ±e ±Si 



±± {s u s 2 ) 



1 + e ±2 ( Sl+S2 ) tanh 2 
(sinh 2si — sinh 2s 2 ) In — 



+4 arctan ( e ±(si+S2) tanh gl - §2 



I _|_ e ±2(si+s 2 ) 

+ 0(y/e), 



(5.30) 



# ±T (si,s 2 ) = +2e ± ( sl - S2 ) + C»(v^), 
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^ ±> (si,s 2 ) = 2 v ^e ±Sl sinhs 2 + 0(e), (5.31) 

e ± <(s 1 ,s 2 ) = zpA-iJL + ^i), 

y/e coshs 2 
6 >± (s 1 ,s 2 ) = -2 v ^e ±S2 sinhs 1 + C(e), 

e<±( 8l ,8 2 ) = ±A_^_ + (i). 

Ve cosnsi 

It is important to note that, in what concerns 6 ±:iz , the equation (5.30) is not the correct 
answer, since in this region all terms Xn i n (5-24) contribute to the 0(1) order. This was 
pointed out by Maldacena and Swanson [28], who derived an integral representation for 
the scattering phase in this region. 

It is interesting to note that ±=fc and 9 ±T in eq. (5.30) depend naturally on the 
rotated coordinates 

S = S\ + S 2 , S = S\ — S 2 , 

and that for large values of s, positive or negative, we have 

# ++ (si, s 2 ) ^ 2e s (tanh ^ - 1 - sinh s) . 

6 Solving the integral equations in the strong cou- 
pling limit with the dressing kernel 

Once we have obtained the leading term of the dressed kernels, we are able to compute the 
leading contribution to the energy of the states we are interested in. In the sw(l|l) and 
su(2) sectors, we do not need to obtain the functional form for the density explicitly. 
This is a fortunate situation, since in the (±±) regime, where most of the roots are 
concentrated, we do not know the kernel in a closed form. The computation is very similar 
to the one performed in [11], where the A 1 / 4 behavior for the anomalous dimensions was 
first derived. 



6.1 The sw(l|l) dressed case 

The results of the numerical and analytical computation by Beccaria et al. [39, 40] for 
the energy of the highest excited state in the su(l|l) sector can be also derived from the 
integral equation. As it was shown in [40], the same results can be obtained using the 
light-cone Bethe ansatz equations [55], first derived in the s«(l|l) sector by Arutynov 
and Frolov [56]. In this sector, the total kernel is 

}C = K + K d . (6.1) 

First, we can show that the densities p< and p> are subdominant, while the main con- 
tribution comes from p±. The density in the "giant magnons" region is given by the 
equation 

P < = ^ + K<< P" + ^ V + E ^ V • (6-2) 
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By symmetry, we have p+(a) = p~{— a) and since 

poo pO 

/ ap + (a)da + / ap^{a)da = , 

Jo J~oo 

we conclude that the ± terms in (6.2) cancel 

P< = ^| + /C<V + /C<V + 0(1). (6.3) 

All the terms except IC << p < are of the order 1. Since /C << = 0(l/e) we deduce that 
p<(s) = 0(e). Similarly, we can write 



p > = + /C><p< + /C»p> + V /C> V . (6.4) 

Z7T as 



Again the leading term in the sum vanishes and all the terms in the r.h.s. are of the order 
e, so that p > (s) is also of order at most e. Therefore the roots of the Bethe equations 
must concentrate in the regions ±, or "near-flat space" region. In this region, we have 

p± = J-!?? + /C±<p< + K, ±> p > (s') + K ±± p ± + /C^pT . (6 .5) 
2tt da 

At the leading order, we are left with 

p ± (a) = J da'K- ±± (a, a')p ± {a) + J da'JC ±T (a,a')p T (a') (6.6) 
or, by specializing to the sign + and by integrating over a, 

oo 

fc+( a ) = __L J da' 9++ {a, a') p + (a') - J da' 9^(a,a')p-(a') , (6.7) 



'tot 



where k + (a) is the counting function, k + (a) = J Q a dap + (a). 9^ Q f is defined as ^ 
— 2n f Q a dct /C +± and at leading order is equal to the dressing part 9 +± . By multiplication 
with p(a) and integration we get 

poo -i poo poo 

/ dak + {a)p + {a) = / da da'9 ++ (a,a')p + (a)p + (a') (6.8) 

Jo 271 Jo Jo 

poo pO 

— — / da da'9 + ~(a,a')p + (a)p~(a') . 
2tt Jo > /_ 00 

Here we used that 9 + ~(0,a') = 9 ++ (0,a') = 0. Let us remind that a = corresponds 
to u — oo. Due to the anti-symmetry of the phase 9 ++ , the first integral in the r.h.s. of 
(6.8) vanishes, while the l.h.s. is equal to 

/•oo pl/2 i 

/ da k + {a)p + {a) = dk + k + = - . (6.9) 
Jo Jo ° 
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Finally, we use that 6 + (a, a') = 2aa' and the symmetry p + (a) = p (—a) to cast the 
equation (6.8) in the form 

J da a p + (a) = • (6.10) 
This is all we need to compute the leading order in the energy 

//»oo 
e-P' 2 sinp/2 p(p) dp ~ Ay/gL J da a p + (a) = y^hTgL . (6.11) 

It is remarkable that the energy of this state is of the order A 1 / 4 . It would be interesting 
to know if there are states in the sw(l|l) sector which have greater energy. This would be 
possible for states where the fraction of giant magnons is larger than 1/ yfg. Presumably 
the strong repulsion among the giant magnons will prevent this phenomenon to happen. 

6.2 The su(2) dressed case 

The computation of the leading order in the energy for su(2) sector goes similarly to 
that for su(l|l), with two differences. First, there is a non-trivial contribution from the 
su(2) kernel to the leading order of 0££ 

u 

f f°° 2e , 

Q su (2)(u, u') = 2tt J dv K su{2 )(v, u') = 2 J dv _ - — = it + 6 asym (u, u) . 

oo 

(6.12) 

Therefore, the first term in the r.h.s. of the equation (6.8) does not cancel anymore and 
we obtain 

J da k + {a)p + {a) = — - f J da p + {a)\ ~\ — I J daap + {a)\ . (6.13) 

Second, the maximum of counting function is 1/4 instead of 1/2, because the total 
number of magnons is M = L/2, and 



PCX) 

/ da p + (a) 
Jo 

From equation (6.13) we obtain 



1/4 1 
dk + = ~ 

4 



J da a p + (a) = (6.14) 



which reproduces 15 the value of the energy obtained by Beccaria et al. [39, 40] 

J su{2) 



E d su(2) = ^L + 0{g»). (6.15) 



15 We thank M. Beccaria for an email exchange which helped us to remedy a discrepancy with their 
results in the first version of this paper. 
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6.3 The si (2) dressed case 

For the sl(2) case, the total kernel with the su(2) piece subtracted, as defined in (2.11), 
is equal to 

K{u, u) = 2K(u, u) + K d (u, u) = —d u <p(u, u) . (6.16) 

The BES equation for this sector takes the simplest form when written for the shifted 
density function 

a(u) = a(u)- — , (6.17) 

7T 

which, up to a negative factor, equals the total density at the leading order in S, see 
(2.16). The BES equation reads 

a(u) = i jH du' j- 4g2 + J du' £(«,«') a(u') . (6.18) 

We are interested in expanding equation (6.18) in powers of e. At small e 

— -T du' ( - a<KU '\ 9 + 0(e 2 )] + [ du' K(u,u') cr(u') = . (6.19) 

* J-oo \[U-uy J J 

This equation splits into several equations, which couple the different regions in u. For 
instance, if \u\ < 1 we have 

I ( £ M 0k + L du ' 0k + o(e 2) ) + L du ' K<< ( "° s< ( "' ) 

+ y rfw' /Co (u, u')a> (u) + J da /C <± (m, a) ^(a) = . (6.20) 

Equation (6.18) was analyzed in the first two orders in e by the authors of [25]. 
They did not consider the regions around u — ±1, which can contain a fraction of the 
roots. Equation (6.18) suggests that the e-expansion of the density starts at order e. 
The structure of the leading terms (5.28) in the expansion of the dressing kernel, where 
half-integer powers of e appear, suggests that the density and the anomalous dimension 
may involve a half-integer powers of e as well. However, these corrections seem to appear 
in higher orders in the anomalous dimension, and this may explain why they have not 
been seen yet. 16 In the following, we are going to ignore the corrections coming from the 
vicinity of the points u — ±1 and solve perturbatively the equation (6.18), by considering 
the decomposition of the density and the kernel in integer powers of e 

a(u) = a(u) = eoi(u) + t 2 <J 2 (u) + e 3 a 3 (u) + . . . , 

7T 

K(u,u') = e- 1 JC_ 1 (u,u') + }C (u,u') + e}C 1 (u,u') + ... . (6.21) 



16 As far we can see, the first non-zero order for a ± (a) is at most e 2 , and there may be corrections of 
order e 5 / 2 to a < (u) and of order e 3//2 for a > (u). The last of these corrections would induce a term of 
order g 1 / 2 in the anomalous dimension. To decide whether they are here or not we have to go to higher 
orders in the expansion of the kernel. 
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• order e° 

At the leading order in e equation (6.20) reads, in simplified notations, 

tC-i&i = 0. 

To this order only the << sector contributes 



1 

du' Kt< (u, u') a$ W) = . (6.22) 

i 



It was proven in [25] that the kernel K,t\ is non-degenerate, therefore we have <jf(u) = 0, 
or equivalently o\ = 2/ir. 

• order e 1 

The 0(e) term in the expansion of equation (6.20) relates o\ and a 2 

/C<< a< +lC<>a> = . (6.23) 
The kernel in the second term is antisymmetric in u' — > —u' 



u 



2 



K%>(u,u') = —d u ^== sgnu' , (6.24) 

7T y/u' 2 — 1 

so that it vanishes upon integration against the symmetric function erf (u). As pointed 
out in [25], the subleading density (u) vanishes as well. It is worth noticing that 
JC^lujii') is the sum of two terms, one coming from the s«(l|l) kernel and the other 
from the dressing phase, /Qp (u,u') = 2Kq > {u, u') + K^q (u,u'). 

• order e 2 

At this order we have an equation relating a±, a 2 and a 3 



2 r 



du'l — — — + £-l°3 + ^0^2 + /CifTi = . 

(u — uy 



o The interval \u\ > 1 

Consider first the region \u\ > 1, where the equation takes the form 

- -f du' ^^'l- + [ du' K,\>(u,u') <J>{u') =0, (6.25) 

71 J > (it - U'Y J > 

with the kernel ^^(u, u') = ^d u (pi given by (5.29 ). Equation (6.25) is a total deriva- 
tive, so we first integrate it to 



£ 



du'a>(u') (---I— + JLfc^u')] =0. (6.26) 

y 71 U — U 2,71 J 
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As the phase is anti-symmetric, there is no integration constant. 

At this point, we find it more useful to switch to the variable s, defined by u = coth s, 
and (6.26) becomes 

1 f°° s — s' 

- 4- ds' a\ (s') sinh s' coth = . (6.27) 

7T J-oo 2 

Now we can reformulate the integral equation (6.26) as a Riemann boundary value 
problem. We first rewrite (6.27) in terms of the normalizable density o\ = o\ + 2/tt 

1 r 00 s — s' 4s 

- + ds' erf (s') sinh s' coth = — — . (6.28) 

7T iy J 2 7T 2 sinh s v ; 

Then, after the redefinition 

Lu = e s , r(u) = -r(l/u) = <T?(s)swhs, (6.29) 
the l.h.s. of (6.28) takes the form of a Cauchy integral 

f^ rM+2 f i/K^US J^_, (u>0) . (6.30) 

The first term on the l.h.s. is actually zero due to the anti-symmetry of r{uj). 

Equation (6.30) can be formulated in terms of a boundary condition for the resolvent 

R(u) = R(l/u) = / dJ -L^lL , (6.31) 
Jo UJ — u 

which has a cut on the positive axis, namely, 

R(u + iO) + R{u - i0) = - lna; _ 1 (cu>0). (6.32) 
The most general solution of (6.32) with the symmetry R(uj) = R{l/u) is of the form 

= 4 M-,) + - V=i7^Q(iM 

7T — 1/UJ 

where Q(w) is a rational function. The latter is determined by the analytical properties 
of the resolvent (no poles outside the positive real axis) and the requirement that the 
density U\ is normalizable 

f f°° f°° 2r(cu) 

/ duaAu) — / dsa^(s) — / du — = finite. 

J> J-oo Jo U ~ 1 

The only solution with these properties corresponds to Q = —1/2: 

4 ln(-cj) - \ (^/^ - yJ-l/J) 
R(u) = ^ ; J - , (6.34) 
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and its discontinuity along the positive real axis, 

r M = ^- = ~ 7. 1 /. . " — -r > ( 6 - 35 ) 



27T2 7T \C<J — C<J — 1 / ' 

reproduces the solution found by Alday et al [25]: 

a< (u) = eof = a > (u) = ea?(u) = (l - cosh . (6.36) 

The total integral of the density, 

r°° 1 7r-4 1 

/ ds(a<(s) + a>(s)) = — + - = -, (6.37) 

J-oo *g Ang 4g 

reproduces the correct leading behavior for the twist-two anomalous dimension 

f(g) = 4g + ... . (6.38) 



o The interval \u\ < 1 

Let us assume for the moment that the corrections from the vicinity of the points u — ±1 
can be neglected. Then, in the region \u\ < 1, the equation will reduce to 



- / du' / 1>(M 1 + K<<a< + K<>&> = 0. 
tJi«'|>i [u-u'y 



(6.39) 



VI: 

The kernels in the first and the third term combine to 



.,2 1 1 „ sinhs'coshs + -is - s') sinh 2 s' 
JC<>(u,u') + — j^ = -d u f± . 6.40 

7T [U — U') 2 7T COsh(s — S f ) 

The procedure to obtain the previous expression was the following: we expressed the 
first two terms xo an d Xi m the dressing phase in terms of the variables p and f3' and 
expanded them up to the desired order in e. Then we expressed the results in terms of 
the variables s, s' by using the relations 

u = cosp/2 = tanhs , v! = cosh/572 = coths' . 

Integrating the product of the kernel (6.40) with erf = — 2cosh(s/2) gives zero, so we 
deduce that 

<7 3 <(«)=0. (6.41) 

• order e 3 

At this order, equation (6.20) gives 



2 f 

7T J-< 



du' G2 U + /C_i<t 4 + JC a 3 + /Ci<7 2 + /C 2 (7i = . 
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The subleading density a\ (u), which we need in order to compute the subleading term 
in /(<?), is supposed to be determined by the equation with \u\ > 1 

- / du ' / 2>(m 1 + k>i > ^ > = -/cr • (6.42) 

After integrating with respect to u and passing to the s parametrization, this equation 
takes the same form as (6.28), but with different inhomogeneous term 

i r°° s — s' i r°° 

- f ds'a> (s') sinh s' coth = — — / ds'af (s') 0>>(s, /) . (6.43) 

7T J-oo 2 27T sinh S ./-oo 

The inhomogeneous term in the r.h.s. can be computed from the expression 



, 4(sinhssinhs') 2 (s — s') cosh(s — s') — sinh(s — s') 

S ' S ) = ~ r T27 A ' ( 6 - 44 ) 



7r sinh (s — s') 

which, after being integrated with o"f(s) = — | cosh(s/2)|u'(s)| gives 
1 



27T sinh s 



/OO 9 
ds' a> (s') 0>> (s, s') = sinh s sinh(s/2) . (6.45) 
■oo TT 



After the redefinition (6.29) we again reduce the integral equation to a Riemann boundary 
value problem 



R(u + ,0) + R(u - ,0) = 2 f = . (6 . 46) 

The solution is 



s sinh s cosh(s/2) s sinh s sinh(s/2) 

R[S) = 2^ ' r(s) = 2^ ' 

which gives for the e 2 -correction to the density 

> ssinh(s/2) 
a * (S) = 2^ • 



(6.47) 



This correction to the density grows exponentially at s — > oo and is not normalizable 17 . 
The significance of the result (6.47) is that the higher order corrections to the density 
become large in the vicinity of the points u = ±1, which correspond to s = ±oo. 
Therefore, in order to determine the next correction to the scaling function f(g) we 
have to take into account the contribution coming from the points u — ±1 (the regimes 
±). The equations starting with (6.39) should be corrected accordingly. We leave this 
problem for future work. 



17 Strictly speaking, the maximum value of s to which we should integrate is s max = K/2 ~ — lne. If 
we stop the integral at s ma x, we will obtain a contribution to the universal scaling function of the order 
lne/\/e, which violates the hypothesis that the expansion is in integer power of e. Also, it signals that 
the limit e — > is not uniform, and this is in accordance with the fact that in the near-flat space regime 
all the terms \n contribute [28]. 
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7 Outlook 



We presented a method to treat the strong coupling limit of the Bethe ansatz equations, 
which can be used to extract the spectrum of anomalous dimensions. We found it useful 
to express the spectral parameter appearing in the Bethe ansatz equations by means of 
a complexified momentum, with real part p and imaginary part f3 

x ± = e m±i P /\ (7.i) 

It is worth noticing that these variables already appeared in [37], where the BDS magnons 
were interpreted as bound states of more fundamental, fermionic excitations. There, p is 
the momentum of the bound states, while (5 is related to the "size" of the bound state. 

Furthermore, the functions of p and (5 are naturally expressed in terms of elliptic 
functions. Our elliptic parametrization is a version of that in [10]. In the strong coupling 
limit, the elliptic parametrization degenerates into a hyperbolic one. In order to take 
into account properly the different regimes, we used different expansions of the elliptic 
functions in terms of hyperbolic functions, depending on the value of the the elliptic 
parameter s: around s = 0, K and ±K/2. These regimes correspond to the "plane- 
wave" [32], the "giant magnons" [33, 34] and the "near-flat space" regions, recently 
characterized in [28] , respectively. It is interesting to note that the momentum p and the 
energy e = E/Ag of the excitations have, in the plane wave and near-flat space regime, 
relativistic-like expressions 

p = esinhs e = e(coshs — 1) (plane wave) (7.2) 
p = ±2^1 e ±s e = 2^1 e ±s (near flat space) (7.3) 

with s playing the role of the rapidity variable. However, only the positive branch of 
the energy, corresponding to particles, appear in the sectors we consider, sectors which 
should be closed at any order in perturbation theory. It is important to understand 
how the antiparticle branch will appear; most probably by another copy of the sector 
which joins in. An example of how this should happen is offered by the su(2) principal 
chiral model [57, 58], where the particles and antiparticles correspond to the two copies 
of sii(2)-symmetric excitations. 

The strong coupling limit of the kernels without the dressing factor is relatively 
simple. The giant magnons interact through a delta term, typical for the statistical 
repulsion. In the sw(l|l) sector, there is also a "length-changing" term. The excitations 
of the plane-wave type do not interact with each other, at leading order, but they interact 
with the giant magnons. The excitations in the near-flat space regime do not contribute 
to the leading term. 

The strong coupling limit of the dressing kernel is more involved. For the giant 
magnons and plane-wave regions, we can do it safely with the strong coupling expansion 
of the dressing kernel [11, 12, 13], at least for the first few orders in e. For the near-flat 
space limit, however, as it was pointed out in [28], all the terms in the strong coupling 
series contribute to the leading order. 

When we consider the kernels with the dressing phase, the repulsion of the giant 
magnons is so large that their density vanishes in the leading order. For sw(l|l) and su(2) 
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sectors, most of the magnons are concentrated in the near-flat space regime. Fortunately, 
although we do not have an explicit expression of the kernel in this regime, it is possible 
to obtain the energy for the highest excited state at leading order. 

The most exciting application of our method is to obtain the strong coupling ex- 
pansion of the twist-two operator anomalous dimension and to compare with the string 
computations [21, 22]. Several recent analytical and numerical works [23, 24, 25] were 
devoted to this task. In particular, [25] obtained the leading term, using the Fourier 
transform representation. They also tried to obtain the result by a method very similar 
to ours, and found that the first two orders of the equation are not sufficient to fix the 
leading order. Here, we show that it is possible to obtained the leading order, by going 
an order higher. We also show that, in order to compute the higher correction, we have 
to take into account the contribution from the near-flat space regime. 

Note added. An integral representation of the dressing phase, seemingly related to 
(1.5), was used in [59] to explore the analytic properties of the S'-matrix . 
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A The elliptic parametrization 

The modulus of the elliptic map is 



The real and imaginary parts of the complexified momentum 

p( s ) = 2k' P(s)=\n±±4i-2 csvdv (A.2) 
Jo dn v 1-k' J K 

have the symmetries 

p(2K-u) = 2n-p(u), P{2K - u) = P{u) 

p(u) = -p(-u), P(-u) = P(u) =F 2in . (A.3) 



The phase in the Bethe equations depends on p and p through 

1 i 

cosh \P = — ns s, sinh \3 = — ds s, 
k z k 
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cos \p = cd s, sin \p = k' sd s. (A.4) 

e ±/j/ 2 = l±dns e ±i P /2 = cds± ik , sdg 
ksns 

d s p = 2k'nds, d s (3 =-2 ess. (A.5) 
Below we give asymptotic expressions of these functions in the limit k' ~ e — > 0, where 

/I 2 

X « In -, K' q = e - nK/K ' « — . (A.6) 

e 2' y 16 v ; 

In this limit it is appropriate to use the expansions of the three elliptic functions in terms 
of hyperbolic functions of v — ns/2K': 

tt / ^ q 2n 

sns = — — — - tanhf + 4> (— ) n — sinh2nv 

2kK' \ ^ V ; 1 + q 2n 



ens = 

2kK 



(1 00 ~2n— 1 N 

_^ + 4V(-) n - 9 - - cosh(2n - 1> 
coshw ^ v 7 1 + g 2 ™- 1 v 7 

/ i 00 ~2n— 1 \ 

dns = —( — 4 V(-) n 9 - - cosh(2n- l)v . (A.7) 

2^'lcosht; ^ v ; l-o 2 ™- 1 v M v ; 



tt / 1 . a 2 "" 1 

x n=l 

In the limit e — > and up to o(e 4 ) terms the three functions are given by [45] 



sns ~ tanh s + e 



2 (sinh s cosh s — s) 



4 cosh s 



1 2 tanh s I 
ens ~ : — — e sinh s 



cosh s 4 V cosh s 

1 9 tanhs / . , s \ , . 

dn s ~ — + e z — - — smh s H — . (A. 8) 

cosh s 4 V cosh s / 

We used only the leading term; it gives a good approximation for s G /> = [— 
Similarly, for the interval s G /< we can use the approximation 

sn (K — s) ~ cds ~ 1 

cn (K — s) ~ fc'sd s ~ e sinh s 

dn (K — s) ~ fc'nd s ~ e cosh s . (A. 9) 

For the functions that enter in the definition of the kernel we find: 
(a) In the interval I> : 

cosh(/3/2) = coths, sinh(/3/2) = P' 



sinh s ' sinh s ' 

cos(p/2) = 1, sin(p/2) = k' sinh s, d s p = 2k' cosh, s 
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1 

u = s - , u = coths 

sinh s 



e^ 2 = coth|, e° = ^^ (A.10) 

(b) In the interval /<, after replacing s — > K — s: 

cosh(/3/2) = 1, sin(/3/2) = fc' cosh s, <9 s/ 3 = -2 A;' sinh s , 

cos(p/2) = tanh s, sin(p/2) = 1/coshs, <9 s p = 2/coshs 



1 

cosh^ s 



u = , u = tanh s, 



Plugging these expressions in (3.27) we evaluate the kernel in the four possible regimes. 
In the sector s, Si G /> the numerator in (3.27) is of order e, while the denominator 
remains finite. Therefore 

K»(s, Sl ) = 0. (A.12) 
For the non-diagonal elements of (3.29) we obtain 

K><(s,s 1 ) = K(s,K-s 1 ) /3'sin(p 1 /2) 



4tt cosh(/3/2) - cos(pi/2) 



27r cosh(s — Si) ' 

T , , x ^/ T , \ 1,1 p' sinh(/3i/2) 
iir<>(s,Si) = iir(K-s,si) ~— v 



4vr 4vr cosh(/3i/2) - cos(p/2) 

g-L^— L / , (A.13) 
27r coshs 27r cosh(s — Si) 

Finally, if both arguments are in I<, then the kernel K«(s, Si) = fT(K — s, K — si) 
vanishes except near the double pole of the denominator at s = s±, where it can be 
approximated by a delta-function: 

K«(s,s 1 )=K(K-s,K-s 1 ) -—p'- 1 p,(3 



4?r 2vr /3 2 + sin 2 (p-pi)/2 
^-p - \d s p\8 (p-pi) 
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2% cosh s 



- 5(s - si). 



(A.14) 



• The regime s ~ ±K/2 

Assume that s = ±K/2 + y, where y <C y/e. We first evaluate the three basic Jacobi 
elliptic functions for the shifted argument: 

sn (s ± ±K) fc'sd ± cn 

sn |K cn 2 + fc'sn 2 

cn (s ± |K) cn =F sn dn ^ 

cn |K cn 2 + fc'sn 2 

dn ( g± l K ) dn ^(l-fcOsncn 
dn±K ~ cn 2 + A;'sn 2 ^' 

When fc' e — > 0, the values of the three functions at s = are 



sn(iK) = 1/VTTk' ~ 1 - |e, 
cn (±K) = v/fc'/(l + fc') « Ve, 



dn(±K) = VF ^v^, (A.16) 



so that 



sn(±±K + y) = ±(l-| ee ^) 
cn(±±K + y) = e TS/ (1 -esinh 2 y) 

dn(±|K + y) = yft e Ty (1 ± e sinh y cosh (A. 17) 

From here one finds for the asymptotics of the functions p(s), /5(s) and w(s) 
for s = ±K/2 + y: 

sinh(/3/2) = ±v^e TJ/ (1 + |e cosh2y) 
sin(p/2) = ±v^e ±J/ (l - \t cosh2y) 

u = ±1 -esinh2y + o(e 2 ). (A.18) 

B Evaluation of K% 

Here we will evaluate the integral (5.21) for the strong coupling limit of K d . We have 

K p _d Pl + K^dfr 



4tt 



K-(xi,x) = 

4:71 

K+(X , X2) = KiE±J^ 
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dz = \Z p dp+\Z p d(3, 



where 

cos(p-pi) -e-l^ 1 ! * _ sin fa - p) 

sinh^cos™ a sin 2^ cosh « 
= ^7 AT£ = ^-7 

Z p = sinh(/3/2)cos(p/2), Z' 3 = sin(p/2) cosh(/3/2) ; 

# a = cosh(/3 + /3 a ) -cos(p-p a ), a = 1,2. (B.2) 

The integral to evaluate is 

i p'Ity poo 

K° d ( Xl ,x 2 ) = -^rj dp J dPiKLdpi + KtdfhHZrKP-ZfiK*) 

1 p2n poo 

= 4^ I dp J d Pi A (P:P)dpi + B( P: P)dp 1 ], (B.3) 

where 

A(J3,p) = (cos(p- Pl )-e-\P + M)G((3,p), 

B((3,p) = -sin(pi -p) G((3,p), 

sin( P -f)sinh(/j+f) + sinhf sin f 
KP,P) " (cosh(/3 + /3 1 )-cos(p 1 -p))(cosh(/3 + /3 2 )-cos(p-p 2 )) ' K ' } 

The strategy of calculation is the following: first we take an integral over p, which is 
of the form: 

I — J dpR(cosp, sinp), (B.5) 
o 

where R is a rational function. By symmetrization p — > — p we reduce the integral to 

J = y dpR(cosp, sin 2 p) (B.6) 
o 

and, after substituting cosp = —t, it can be evaluated by taking the residues. 

The integral over (3 becomes, after the substitution (3 = — logx, an integral from 
to 1 of a rational expression and can also be performed explicitly. As a final result, we 
obtain 

oo 2n , \ 

11 1 1 



M 





dp A(p, f3) = l n\— + -(xt-—) log - c.c. ) , (B.7) 
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OO 2lT / j 



J dp J dpB(p, (3) = 7T 1 ±= + \{xt - ^p) log 1 | + c-c. J ■ (B.8) 





From here we deduce the expression for K®, presented in (5.22). 
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